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Abstract 

When a sinusoidal (ac) field is applied to a suspension containing nonlinear 
dielectric particles, the electrical response will generally consist of ac fields at 
frequencies of the higher-order harmonics. The situation is further compli- 
cated by an intrinsic dielectric dispersion which often occurs due to the surface 
conductivity or inhomogeneous structure of the particles. We perform a per- 
turbation method to investigate the effect of intrinsic dielectric dispersion on 
the harmonics of local field as well as induced dipole moment. The results 
showed, for weak intrinsic dispersion strength, the ratio of the third to first 
harmonics of the induced dipole moment decreases as the frequency increases, 
which is qualitatively in agreement with experimental result. However, for a 
strong dispersion strength, the harmonics ratio increases as the frequency in- 
creases. Moreover, an increase in the intrinsic relaxation time may increase 
the strength of harmonics. 

PACS Number(s): 72.20.Ht, 83.80.Gv, 82.70.Dd, 41.20.-q 
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I. INTRODUCTION 



When an intense electric field is applied to a suspension containing highly polarized 
particled embedded in a host fluid, the particles will form chains due to the induced dipole 
moments inside the particles, and hence complex anisotropic structures occur. Also, the 
apparent viscosity of the suspension will be enhanced greatly. These structures can have 
anisotropic physical properties, such as the effective conductivity and permittivity [1,2]. In 
this case, optical nonlinearity enhancement was found as well [3]. 

A convenient method of probing the nonlinear characteristics of the fleld-induced struc- 
tures is to measure the harmonics of the nonlinear polarization under the application of a 
sinusoidal (ac) electric field [4]. Therefore it is possible to perform a real-time monitoring of 
the field-induced aggregation process in such a suspension by measuring the optical nonlin- 
ear ac responses both parallel and perpendicular to the uniaxial anisotropic axis. Recently, 
an anisotropic Maxwell-Garnett theory [5] was applied to study the nonlinear polarization of 
these anisotropic structures [6] . Anisotropic optical nonlinear ac response was indeed found. 

When a nonlinear composite with nonlinear dielectric particles embedded in a host 
medium, or with a nonlinear host medium is subjected to a sinusoidal field, the electri- 
cal response in the composite will in general be a superposition of many sinusoidal functions 
[7]. It is natural to investigate the effects of a nonlinear characteristics on the optical non- 
linear ac response in a colloidal suspension which can be regarded as a nonlinear composite 
medium [8,9]. The strength of the nonlinear polarization should be reflected in the mag- 
nitude of the harmonics. For these field-induced anisotropic structures, we found that as 
more and more chains are formed the harmonics of local field and induced dipole moment 
increase (decrease) for longitudinal (transverse) field case [6]. 

On the other hand, the effect of nonlinear characteristics on the inter-particle force has 
been analyzed in a colloidal suspension of nonlinear particles [10] and further extended to a 
nonlinear host medium [11]. 

Actually, the intrinsic dielectric dispersion often occurs due to the surface conductivity 
or the inhomogeneous structure (e.g., coating effect [12]) which often exist in colloidal sus- 
pensions. Hence it is of particular importance to discuss the effect of an intrinsic dielectric 
dispersion on the nonlinear polarization. In this regard, we will use a perturbation expansion 
method to calculate the nonlinear ac response of a nonlinear composite, i.e., harmonics of 
local field as well as induced dipole moment. 

II. NONLINEAR POLARIZATION AND ITS HIGHER HARMONICS 

Let us consider the system in which the suspended particles of radius a with a nonlinear 
dielectric constant are suspended in a host medium of linear dielectric constant 62. The 
nonlinear characteristic gives rise to a field-dependent dielectric coefficient [13]. In this case, 
the electric displacement-electric field relation inside the spheres is given by: 

Di = eiEi + XilEil^Ei ^ eiEi + xi(|Ei]2)Ei 

= {eiH + , )Ei + Xi(|Ei|2)Ei = eiEi, (1) 

1 + lUj/UJo 



2 



with i — V— Tj where ei and Xi the hnear coefficient and the weak nonhnear coefficient 
of the suspended particles, respectively. Here ei is expressed as a Debye expression, 
is the dielectric constant at high frequency, Aei is the intrinsic dispersion of the particle, 
cuq — 2'k/to is the intrinsic characteristic frequency, with tq being the intrinsic relaxation 
time, and cu is the angular frequency of the external field. In Eq.(l), we have adopted an 
approximation: the local field inside the particles is assumed to be uniform. This assumption 
is called the decoupling approximation [14]. It has been shown that such an approximation 
yields a lower bound for the accurate result for the local field [14]. We further assumed 
that xi is independent of frequency, which is a valid assumption for low-frequency processes 
in colloidal suspensions. As a result, the induced dipole moment under an applied field 
E(i) = E{t)z along axis, with t being time, is given by [6]: 

p = iea^bEit), (2) 

where eg is the effective dielectric constant of the suspension, and b is the field-dependent 
dipole factor: 



ei - £2 
ei + 262 





|Ei 


n-c2 


ei + Xi(|Ei| 


2) +262 



(3) 



In order to obtain the effective dielectric constant ie, we invoke the Maxwell-Garnett 
approximation (MGA) for anisotropic composites [5,6]. For the longitudinal field case when 
the ac field is applied along the uniaxial anisotropic axis, the MGA has the form 

— ^2 _ ^ ei — 62 



Pl^s + (3 - Pl)€2 ei + 2e2 

with / being the volume fraction of particles, whereas for an transverse field case when the 
ac field is applied perpendicular to the uniaxial anisotropic axis, the MGA expression is 
obtained by replacing Pl with — (3 — /3l)/2 [5,6]. Here Pl and Pt denote the local field 
factors parallel and perpendicular to the uniaxial anisotropic axis. These local field factors 
are defined as the ratio of the local field in the particles to the Lorentz cavity field [5]. For 
isotropic composites. Pi, = Pt = I, while both /3l and Pt will deviate from unity for an 
anisotropic distribution of particles in composites. These P factors have been evaluated in 
a tetragonal lattice of dipole moments [5] and in various field-structured composites [15]. 

We resort to the spectral representation approach [16] to represent the local field inside 
the particle [3]: 



|Ei| 



f 



r™d., (5) 

Jo \s — xr 



with the material parameter s = (1 — 61/62)^^ as well as the spectral density function 
IJ,{x) = fS[x — (1 — //9)/3], which satisfies the sum rule Jq n{x)dx — f. Hence Eq.(5) has 
the form 
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Actually, this is a self-consistent equation for (|Eip), which can be solved at least numerically 
[6]. In what follows, we will use a perturbation method instead to obtain analytic expressions 
for harmonics of the local electric field and the induced dipole moment. 

On the other hand, due to the uniform local field in each particle, Eq.(6) can also be 
directly derived from [13,17]: 



dei 



(7) 



If we apply a sinusoidal electric field, i.e. E(t) = Esin(a;t), the induced dipole mo- 
ment p will depend on time sinusoidally, too. By virtue of the inversion symmetry, p is a 
superposition of odd-order harmonics such that 



P — Poj sin out + sin 3out 
Also, the local electric field contains similar harmonics 

El 1 2) ^ sin cut + sin ?^t -\- 



(8) 



(9) 



These harmonic coefficients can be extracted from the time dependence of the solution of p 
and Ei{t). 



III. ANALYTIC SOLUTIONS 



In this section, we will apply the perturbation expansion [8,18] to extract the harmonics 
of the local electric field and the induced dipole moment. It is known that the perturbation 
expansion method is applicable to weak nonlinearity only, limited by the convergence of the 
series expansion. 



We expand p and yxidEip) into a Taylor expansion, taking xi(|Eip) as a perturbative 
quantity: 



~ a^Eit) a" 



■ r gi(l + 3/-//3) + e2(2-3/ + //3) 

6i(l-//3) + 62(2 + //?) 



(Xl(|El|^))^ (10) 



ei=ei 



Xi(|EiP) = E 



1 9"^ 



3|e2U/xiE2(t) 



:i-//3)ei + (2 + //3)e2| 



2 \ \ m 



(Xi(lEilO) 



(11) 



In view of weak nonlinearity, we can rewrite Eqs.(lO) and (11), keeping the lowest orders of 
Xi(|Ei|2) and xiE2(i): 



P^hiE{t) + h3XiE%t), 



VxidEf I) ^ ji^xim)+h Uxim)) , 



(12) 



where 



4 



hi — a €2 



3 ei - €2 p + 3/(ei - 62) 



ei + 2e2 



Ji = 



27a3|e2|2ei 


"3/(61-62) 


p|p|2(ei + 2e2) 




3|e2| . 27 
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+ 



\p[ 



61 + 2e2 
1 . 



withp=(l-//?)ei + (2 + //5)e2. 

In the case of a sinusoidal field, we can expand E^{t) in terms of the first and the third 
harmonics. The comparison with Eq.(8) yields the harmonics of the induced dipole moment: 



Similarly, we find the harmonics of the local electric field 



X 



(13) 



(14) 



In the above analysis, we have used the identity sin^a;i = (3/4) sina;^ — (1/4) sin3a;t to 
obtain the first and the third harmonics. Similar analysis can be used to extract the higher- 
order harmonics, by retaining more terms in the series expansion. For instance, to obtain 
fifth harmonics, we may take into account the identity sin^ ut = (5/8) sin cjt— (5/16) sin3u;t+ 
(1/16) sin 5a;t, and then obtain the harmonics of induced dipole moment and local electric 
field, respectively: 

3 5 

= hiE + -h^XiE^ + fi5xlE\ 

1 5 

P3u; = --rhxiE^ - T^KxlE\ 
4 Id 



where 



Sla^leal 



9/e: 



p'iei + 2e2 



p2(l-//3)(6i-62) ^ 



369 



+ 



36i 



61 + 262 P(6l + 262)^ 



(1 + 3/ - fP)p - (1 - //?)(p + 3/61 - 3/62) _ 2p+ 6/(61 -62) .-, 

^ P 61 + 262 



and 



1/2 77. 
Xi E^ 



X^E^ 
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-{kixTEf - ^^j.ix'.^Ef, 



(15) 
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where 

,3 1 1 3 p2 



243(1 -//?)2| 
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4 
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' '"'^ ^2p2|p| + |p|3 ' 2 |p|5'- 

Because of the weak nonhnearity under consideration, the effect of high order harmonics is 
quite smalL Thus, without loss of generahty, we only consider the Taylor expansion up to 
second term. In this regard, we will use Eqs.(13) and (14) to investigate the nonlinear ac 
response in the numerical calculations in next section. 



IV. NUMERICAL RESULTS 

We are now in a position to perform numerical calculations to investigate the dispersion 
effect on the harmonics of local electric field as well as induced dipole moment. We let 
/ = 0.15, EiH — 6eo, £2 = 75eo, where eo is the dielectric constant of free space. Additionally, 
set po = eeO^^ei — e2)/(ei + 2e2)E to normalize the induced dipole moment, where is 
given by the anisotropic MGA (Eq.(4)) without nonlinear characteristic. Obviously, the 
expressions for the normalized local electric field and the normalized induced dipole moment 
are functions of a single variable Xi'^E [9]. In view of the weak nonhnearity, we set Xi'^E — 
0.9y^ for our calculation. 

In Fig.l(a)~(c), we investigate the harmonics of local field and induced dipole moment 
versus frequency for different intrinsic dispersion strengths Aei. It is obvious that increasing 
the dispersion strength reduces the harmonics of both the local field and the induced dipole 
moment. We find that, for increasing frequency, the harmonics of local field always increases. 
However, the harmonics of induced dipole moment exhibits a more complicated behavior. 
For small Aei , the ratio of the third to first harmonics of induced dipole moment decreases 
as the frequency increases, which is qualitatively in good agreement with the experimental 
data [4]. However, for large Aei, its harmonic ratio also increases as the frequency increases. 
We expect this behavior can be demonstrated by experiment in the future. 

In Fig. 2, we investigate the harmonics of local field and induced dipole moment versus 
frequency for different intrinsic relaxation times Tq. A longitudinal field case {(Sl = 2) is 
discussed in Fig.2(a), while a transverse field case {Pt = 0.5) in Fig.2(b). It is shown that, 
in either case, increasing relaxation time tq may enhance the harmonics of both local field 
and induced dipole moment. In addition, the transverse field case may predict stronger 
harmonics than the longitudinal field case within high-frequency region, which is actually 
valid for different dispersion strengths as well (no figures shown herein). 

On the other hand, Klingenberg claimed that increasing external electric field leads to 
increasing harmonics of electric current [4]. In this sense, we have also calculated the effect 
of external electric field on the nonlinear ac response of the induced dipole moment (no fig- 
ures shown herein). Obviously, our calculations qualitatively demonstrate the experimental 
results [4] as well. 



V. DISCUSSION AND CONCLUSION 

Here a few comments on our results are in order. In the present study, we have examined 
the case of nonlinear particles suspending in a linear host. Our considerations may be 
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extended to a nonlinear host medium [11]. 

We have studied the effect of an intrinsic dispersion on the nonhnear ac response of a 
colloidal suspension. Actually, our formalism can be readily generalized to discuss two or 
more intrinsic dielectric dispersions of particles (or host fluid) . In view of the fact that the 
intrinsic dielectric dispersion often occurs due to the surface conductivity or the inhomoge- 
neous structure (e.g., coating effect [12]), we believe such effect also plays an important role 
in the dielectric behavior of biological cells, such as the electrorotation, dielectrophoresis 
and electro-orientation. 

The perturbation method is used in the present work, but we believe similar results may 
also be predicted by a self-consistent theory [6] . 

In summary, a perturbation method has been employed to compute the local electric 
field and the induced dipole moment for suspensions in which the suspended particles have 
a nonlinear characteristic and an intrinsic dispersion, in an attempt to investigate the fre- 
quency effect on the nonlinear ac response. The results showed, for weak intrinsic dispersion 
strength, the ratio of the third to first harmonics of the induced dipole moment decreases 
as the frequency increases, which is qualitatively in agreement with experimental result. 
However, for a strong dispersion strength, the harmonics ratio increases as the frequency 
increases. Moreover, an increase in the intrinsic relaxation time may increase the strength 
of harmonics. 
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FIGURES 



FIG. 1. Harmonics of the local electric field and the induced dipole moment versus frequency 
for To = W~^-^s and Pl = 2, for different intrinsic dispersion strength, (a) Aei = 50eo; (b) 
Aei = lOOeo; (c) Aei = 150eo. 

FIG. 2. Harmonics of the local electric field and the induced dipole moment versus frequency 
for Aei = 50eo, for different intrinsic relaxation time tq. (a) = 2; (b) Pt = 0.5. 
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